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Entanglement is a ubiquitous feature of low temperature systems and believed to be highly 
relevant for the dynamics of condensed matter properties and quantum computation even at 
higher temperatures. The experimental certification of this paradigmatic quantum effect in 
macroscopic high temperature systems is constrained by the limited access to the quantum 
state of the system. In this paper we show how macroscopic observables beyond the energy 
of the system can be exploited as proxy witnesses for entanglement detection. Using linear 
and semi-definite relaxations we show that all previous approaches to this problem can be 
outperformed by our proxies, i.e. entanglement can be certified at higher temperatures 
without access to any local observable. For an efficient computation of proxy witnesses one 
can resort to a generalized grand canonical ensemble, enabling entanglement certification 
even in complex systems with macroscopic particle numbers. 
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I. INTRODUCTION 

While the occurrence and possible uses of entanglement were first studied for bipartite states, 
entanglement in systems containing a large number of particles is of interest both from a theoret¬ 
ical and from a practical point of view. Even though the macroscopic world we experience daily, 
can be described classically, there are a number of systems that are large enough to be described 
by the thermodynamic limit, which exhibit quantum behaviour, Bose-Einstein condensates, fer¬ 
romagnetic and superconducting materials being prominent examples. Entanglement may turn 
out useful in understanding thermodynamic phenomena such as phase transitions in such systems 
m- Recently there has also been a lot of attention on the role of entanglement in quantum ther¬ 
modynamics 11H6]. Other possible applications of large entangled systems are quantum computers 
based on solid state or NMR systems [ZHini. In addition to studying entanglement in the limit 
of many particles it is also worth asking up to which temperature entanglement can exist. This 
is an important question for experiments, where cooling down systems requires lots of resources. 
While entanglement usually exists at very small temperatures, it could persist to up to lOOK in 
superconductors [m. 

Experimentally detecting entanglement in macroscopic systems is generally a highly non-trivial 
task. Checking for instance the famous PPT (positivity under partial transpose) criterion, as easy 
as it is theoretically, requires a full state tomography, which is not possible in large systems. Also, 
calculating the eigenvalues for matrices of large dimensions is not practical. The method of choice 
are entanglement witnesses, i.e. observables with positive expectation value for all separable states 
but with negative expectation value for some entangled states. Witnesses reduce the complexity of 
entanglement detection to the measurement of a single observable. However, this observable might 
have no physical meaning and might be hard or impossible to measure. In particular it might be 
necessary to perform a collective measurement of all particles, which is not experimentally feasible 
in macroscopic systems. What is feasible is the measurement of macroscopic observables such as the 
mean energy, the magnetisation, the temperature or the entropy of the system. There have been 
several results showing that mean energy and temperature can serve as entanglement witnesses at 
low temperatures ([IMTj to name just a few). 

However, all of these are limited intrinsically at higher temperatures when the value of the 
macroscopic witness is consistent with separable pure states. Knowing that any experimental 
system will have some non-zero entropy, i.e. not be in pure state, often allows for reasonable 
lower bounds on the system entropy to be assumed (as for example through ambient temperature 
and the second law of thermodynamics). Here we can leverage the entropy to bound “by proxy” 
generic entanglement criteria, allowing us to detect, in principle, all entangled Gibbs states of a 
many-body system, as well as entangled states far out of equilibrium. This method can in general 
be phrased as a semi-definite program (SDP), which are efficiently solvable for small system sizes 
(and have in fact often been used in the context of entanglement quantification [1814201 1. In the 
following we will showcase some exemplary situations where these SDPs can improve entanglement 
detection for systems of up to thirteen qubits on a regular laptop. Furthermore we show that using 
entanglement witnesses in particular allows one to harness tools from statistical physics and solve 
the problem through introducing a virtual “chemical potential” in a generalized Gibbs ensemble, 
changing the maximum entropy state for a given energy. In these generalized Gibbs ensembles 
entanglement witnesses play the role of additional conserved quantities, making the proxy method 
as accessible for large systems as the computation of Gibbs entropies, which we demonstrate by 
detecting entanglement by proxy in the thermodynamic limit. 

In what follows we will assume the Hamiltonian to be well characterized, which for large systems 
is of course only an approximation. Unfortunately, to characterize entanglement one requires 
reasonably precise knowledge of the observables used to certify it. As we will show later, the 
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method is robust against small perturbations, but if the Hamiltonian is entirely different it can 
of course lead to false positives. Our method is suitably generic and requires knowledge only of 
conserved quantities for use in the generalized Gibbs ensembles. While we generically use the 
average energy to showcase our methods, they could just as well be replaced by other macroscopic 
approaches such as e.g. spin squeezing [2TH23] . 


II. ENTROPY AS ENTANGLEMENT WITNESS 


A. The primal problem 


As mentioned in the introduction, the mean energy of a system can be used to witness entan¬ 
glement in the corresponding quantum state. Namely any state with mean energy less than 

-E^min.sep = min Tr pH (1) 

pSsep 

is entangled. By convexity the minimum is attained in a pure state, sep can be chosen to be the 
set of fully separable states or the set of /c-separable states. For k = 2, genuinely multipartite 
entanglement is detected. If the system is in thermal equilibrium, it is possible to derive an 
analogous criterion for the temperature. 

The goal of this chapter is to find a condition that is able to detect entanglement at higher 
energies than Hmin^sep- The idea is to add additional constraints to Q. For example, a lower 
bound on the (von Neumann) entropy of p: 

•Fmin,sep,S’ = min Tr pH. (2) 

pSsep, S{p)>S' 

If E{p) < F'min.sep.S’ and S{p) > S', p is entangled. S' can be varied between 0 and Ind. Clearly, 
it holds F^min,sep,S’ ^ .Fmin.sep- But the relevant question is whether there are entropies for which 
the strict inequality holds. While Q is minimised by a pure state, the minimum does not have to 
be unique. Hence a mixture of minimisers could also have mean energy Fimin.sepj but at non-zero 
entropy. The exact behaviour depends, of course, on the Hamiltonian. However, it is possible to 
show that if there exists an S', for which the strict inequality holds, it will hold for any larger 
entropy. This follows from 


Lemma 1 Emin, sep, S’ os function of S' is convex. 

Proof Let 0 < Sip < Ind and 0 < p < 1. Let pip be separable states with entropies 5(2 
respectively, that minimise Then 

-Fmin,sep,S’ {pS'i + (1 - p)S' 2 ) < Tl H (ppi (1 -p)p2) = pEm\n,sep,sfS'i) (1 “ p)£^min,sep,S’ (5'2) , 

( 3 ) 

where the inequality is due to the fact that ppi -|- (1 — p)p 2 is a feasible point of Q, which follows 
from the concavity of the entropy. □ 

Since F'min,sep,S’ is convex, it is strictly monotonically increasing as soon as it exceeds Fimin,sep- Let 
us call the smallest entropy where the constraint hits Smin- See also figure [T] 

Since for S' > S mm , Limin,sep,S’ is strictly monotonically increasing, we can obtain the same 
curve, i.e. detect the same entanglement, by computing 


5, 


max,sep,E 


max Sip) 

pSsep, Tr/fp=E 


( 4 ) 
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FIG. 1. Qualitative form of i5min,sep,S’ as a function of S'. The dark grey area under ifmin.sep is detected 
by [I 2 HI 4 ] . The method presented here can detect entanglement in the light grey area, i.e. for S > S'min 
and E < ifmin.sep.S’- In the white area above the Flmin,sep,S’ curve, separability is possible. Because of it’s 
convexity, the same curve can be obtained by maximising the entropy under the constraint that the mean 
energy is less than E and varying E. 


and varying Emin,sep < E < Equivalently we could also demand Tr Hp < E. Note that if we 

remove the separability constraint @ will become the Gibbs state entropy 


5’Gibbs = max S{p), 

p state, T^vHp=E 


(5) 


attained by the Gibbs state. Since the optimisation is of the same form as the one yielding the 
regular Gibbs ensemble we will call the resulting state a separable Gibbs ensemble. In order to 
compare its entropy with the Gibbs state entropy, let us define the entropy gap AS = Soibbs — 
'S'max,sep,E' 


Theorem 2 For mean energy E and corresponding temperature T, any state p with Smax,sep,E < 
S{p) < Scibbs is entangled. In particular the Gibbs state is entangled if AS > 0. 

Note that for any non-zero gap between 5max,sep,E and 'S'cibbs there always exist states different 
from the Gibbs state which we detect, for example there is always a state with entropy AS'max,sep,E + 
(1 - A)5Gibbs, where 0 < A < 1. 

The optimisation in Q is difficult to deal with because it contains the separability constraint. 
The main idea is to relax this constraint using sets of states which remain positive semidefinite 
after the application of positive (yet not completely positive) maps A. A prominent example would 
be the partial transposition. While these are of course supersets of the separable states, it is clear 
that for every entangled state in principle there exists a map A and thus a semidefinite relaxation 
that will still yield optimal results for the constrained optimisation. 

Let us use the fact that every set of A-positive states forms a convex set that can be approxi¬ 
mated by a suitable set of entanglement witnesses Wi with suitable weights > 0 [2Tj- As every 
optimal entanglement witness for the set of A-positive states for partition A can be written as 
A^ ® 11 ^[|V’)(V’|] 0116 can find suitable entanglement witnesses for the system in question. While 
this of course comes at the expense of finding suitable entanglement witnesses for specific systems, 
it also opens the possibility to constrain the entropy beyond just states which are separable under 



5 



FIG. 2. Schematic view of Corollary for one witness W. The constrained entropy 5'max,iv,£; or short Sw 
is attained at the intersection between the (green) hyperplane induced by the witness Tr(pVF) = 0 and the 
(red) hyperplane induced by the mean energy Tr(piJ) = E. The concentric circles represent equi-entropic 
states, with the outside low entropy region (blue) bounded by Sw and the center high entropy region (red) 
bounded by the Gibbs state entropy Scibbs- The detected states lie within the (white) region spanned by 
AS. 

fixed bi-partitions. I.e. it enables us to find also genuine multipartite entanglement or any other 
non-partially separable set, by choosing corresponding witnesses. In order to remain fully general 
we include both an arbitrary set of entanglement witnesses and positive maps in the following 
considerations. 

The relaxed problem is thus given by 

*S'max,sep,E ^ 5'jnax,A^,Wi,i? =maxS'(p) (6) 

s.t. p state, Tr {pH) = E, 

Aa ® ^aIp] - O’ > 0- 

‘S'max,A^,Wi,£; will be an upper bound on S'max,sep,E, it’s tightness depending on the choice of wit¬ 
nesses or maps. Since the relaxation just provides an upper bound for the separable entropy it 
trivially follows that according to Theorem 

Corollary 3 For mean energy E and corresponding temperature T, any state p with Smax,AA,WiF < 
S{p) < Scibbs is entangled. In particular the Gibbs state is entangled if Smax,AA,WiF < Sgibbs- 

While the relaxation weakens the detection criterion, we are now dealing with a convex optimisation 
problem with only linear and semidefinite constraints (since maximising a concave function is of 
course equivalent to minimising a convex one). Such programs can be solved efficiently numerically 
thanks to so-called interior point methods\lf)\. In addition, they have a duality theory which can 
be used to give certified upper bounds on 5'max,A^,iyi,£;) discuss in the next section. One 

might ask why instead of applying Corollary]^ one cannot simply apply the witnesses or positive 
maps involved directly. While this is possible theoretically, let us note again that witnesses or 
positive maps are in general not easily accessible in experiments, while macroscopic variables such 
as the mean energy and entropy are. Since Corollary makes use of witnesses without the need 
to measuring them directly, but instead requires the measurement of mean energy and entropy, we 
call those two quantities Proxy Witnesses. 

Let us conclude this section by noting that apart from the von Neumann entropy, in principle 
any concave function can serve as a proxy witness. An example is the so called linear entropy, 
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also known as impurity of a quantum state, which is defined by Slip) = 1 — Tr/o^ |26j . The linear 
entropy is upper bounded by the von Neumann entropy and can be seen as a measure of mixedness 
of a quantum state, as well. The optimisation of the linear entropy is easier to deal with than of 
the von Neumann entropy, as it is only quadratic in p. 


B. The dual problem. 

For convex optimisation problems, such as[^ it is possible to derive a dual optimisation problem. 
To do so, one defines the Lagrangian, a function incorporating both objective, i.e. the function to 
be maximised, and the constraints. The constraints are added by means of Lagrange multipliers. 
The Lagrange multipliers are referred to as dual variables, whereas the variables of the original 
problem are referred to as primal. Maximisation of the Lagrangian over all admissible primal 
points yields the dual objective. The dual problem is then given by the minimisation of the dual 
objective with respect to the dual variables. It can be shown that any dual feasible point provides 
an upper bound on the original problem, which is referred to as weak duality. For details, please 
refer to [25]. For|^the Lagrangian is given by 

L{p,\,pL,Vi,XQ,XA,K) = S{p) + A(Trp - 1) + ^(TrpiL - E) + Tr pWj 

* (7) 

+ Tr pXo + ^ Tr Aa 0 1]4[/o]Xa,a, 

A,A 

where A, p and the Vi are real Lagrange dual variables, corresponding to the trace, mean energy 
and witness constraints, respectively. As the witness constraint is given by an inequality, we can 
restrict to > 0. Xq and the Xa,a, the Lagrange dual variables corresponding to the positivity 
and positive maps constraints respectively, are positive semidefinite matrices. 

Here it becomes clear that witnesses are easier to deal with numerically than positive maps, as 
they only require scalar variables. The dual objective function reads as follows 

e{X,p,Ui,Xo,XA,A) ■.= maxL{p,X,p,i^i,Xo,XA,A), ( 8 ) 

p 

where the maximisation is over all admissible p, i.e. those for which the right hand side is defined, 
in particular not necessarily obeying the primal constraints in Q. Note that the dual problem is 
also a convex optimisation problem: the constraints are indeed linear and semidefinite, while i is 
convex, owed to the linearity of L in the dual variables. 

Clearly, for every primal feasible p and dual feasible X,p, t'i, A"o and Xa,a it holds 

S{p) < L{p, X,p, ni,Xo,XA,A) < ^iX,p, ni,Xo,XA,A) (9) 

In particular we have what is referred to as weak duality: 

*S'max,AA,H^^,£; < rami{X, p, l^i, Xq, Xa,a), (10) 

where the minimisation is over the dual feasible region, i.e. A, // G M, > 0, Xq > 0, Xa,a > 0. 
The right hand side is referred to as the Lagrange dual problem. Note that any dual feasible bound 
yields an upper bound on Smax.AA.Vl/i,^;- This allows us to obtain analytical upper bounds from 
numerics. One simply has to numerically optimise the dual problem and check if the so obtained 
optimisers are dual feasible. If they are they can be inserted into i yielding the bound. 

Let us now evaluate the dual objective function Q for the von Neumann entropy, which is well- 
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defined on any positive semidefinite matrix. Since the entropy is concave, so is the Lagrangian. 
Hence it is sufficient to find a critical point, that is positive semidefinite, of L as function of p. 
Using 'VpS{p) = — 1 — Inp and Vp Tr pM = M with respect to the trace inner product on Hermitian 
matrices, as well as the fact that Tr (A^ ® 1.^[p]Xa,a) = Tr (g) 1.-a[Xa,a^ one can obtain 

V pL = —U — lnp + All + pH + 'y ^ ViWi + A^o + ^ ^ A^ 0 ll^[A^^yv] i (H) 

i A,A 


which vanishes for 


= exp (A - 1)1 + ^77 + u,W^ + + aA[^A,A] 


A,A 


( 12 ) 


Hence the dual objective is given by 

7(A, p, V,, Ao, Aa,a) = L{p<^"\ A, p, ui, Ao, Xa,a) (13) 

= Tr exp ^- 1)1 + pH+ Y, + Aq + A^i 0 ^a[^a,a]^ - X - pE 

(14) 


A,A 




Tr exp j pH + ^ i/jlUj + Aq + ^ A^ 0 1^[Xa,a] j - A - pE. 
i A,A 


(15) 


We will have to minimise i with respect the dual variables. Since Aq > 0, the minimum will be 
attained at Aq = 0. i can be easily minimised for A. This is due to the fact that the eigenvalues 
of a hermitian matrix cannot decrease if a positive semi-definite matrix is added, which follows 
from Theorem 4.3.1 of m- By convexity the minimum is attained where the derivatives vanishes, 
yielding a function i just of p, t'j, A^^a- 


i{p, Ui,XA,A) = InTr exp ^ pH i^iWi -7 ^ A ;4 0 1^[Aa,a] j - pE. (16) 

By weak duality, it holds 


A,A 


'S'max.A A .W,, 








(17) 


As i is convex, i is convex, as well [25]. If only witnesses are used, i only has scalar variables, 
which results in a greatly enhanced numerical performance compared to the primal problem. 

It is an interesting observation that i is reminiscent of a grand canonical ensemble and 

i = \nZ'+ /3E = lnTve^p{-/3H') + l3E, (18) 

where H' = H — ^ • ViWi + Y^a a ^a ® ^a[A^-4a]) • Instead of particle numbers the constraints 

stem from the specific witnesses or maps used. We will sometimes refer to eq@ as a witness 
canonical ensemble entropy (WCEE). 

In other words, when choosing p = —/3, the operator A)^01 ^[Aa,a] can be seen 
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as some sort of grand canonical ensemble with additional ’’chemical” or rather witness potentials. 
In particular the entropy gap can be lower bounded as follows: 


AS > In- 

Trexp 


Trexp(-/3/I') 

+ E. + Ea,a ® t^[XA,A]) 


(19) 


for all t'j > 0 and Xa,a > 0. 

It is also possible to compute i for the linear entropy. After rewriting the Lagrangian as 


Alin — 1 — A — laE + Tr p | All + fj,H + ViWi + Xq + j ~ Tr , (20) 

A,A 


it is easy to show that the optimiser is given by 


= 


- j A1 + pH + t'jVFj + Xq + Aa (8) 11^[Aa,a] j • 


( 21 ) 


A,A 


Hence 


^lin(A, p, t'i, Xq, Xa,a) ~ 4 I ^iWi + Xq + A^ ® I1^[Aa,a] j +l — X—pE. (22) 

A,A 


Minimisation of £lin yields an upper bound on the maximal linear entropy achievable by separable 
states. Again, this minimisation is easier to deal with as in the von Neumann case, because it is 
only a quadratic function. 


C. On numerics 


Let us now briefly discuss how the primal and dual optimisation problems introduced in the 
preceding sections can be implemented numerically. As mentioned before, Q and ( |17[ ) have concave 
and convex objectives, respectively, as well as linear and semidefinite constraints. Note that for 
our purposes it is sufficient to only compute the dual problems. It can however be instructive 
to also compute the primal problem in order to obtain the optimiser and check if strong duality 
holds. While interior-point methods can in principle solve such problems [25], readily available 
solvers such as Sedumi |28| or SDPT3 |29| can only handle linear and quadratic objectives. This 
is sufficient to solve Q and minimise (22) for the Linear entropy. 

For the von Neumann entropy, however, there is a way to obtain an approximate solution: It 
is a well known fact that the von Neumann entropy of a state is equal to the Shannon entropy of 
its eigenvalues. As for i, note that for a hermitian n x n matrix M with eigenvalues A* it holds 
Trexp(M) = Hence £ is a function of the eigenvalues of the exponent. Both S and I 

are invariant under permutation of the eigenvalues. This allows us to reformulate ([^ as 


5'max,AA,Wi,£ =maxiL(v), (23) 

s.t. V = eig(p), p state, Tr (piL) = E, 

Aa (8) IaW - Tr(pITj) > 0 


where eig(/9) denotes the vector of eigenvalues of p. Note that this is not a semidefinite constraint. 
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There is, however, a trick to include the eigenvalues into a semidefinite programme [30]: While the 
eigenvalues are generally not SDP-representable, the sum of the k largest eigenvalues of a matrix is 
m- We can now replace the eigenvalue constraint in ( |23[ ) by the constraint that v has to majorise 
eig(p), i.e. 


Tr(p) = 

i 

Slip) < 

S2ip) <Vi+V2 


n—1 

Sn-lip) < 

i=l 

where Sk{M) denotes the sum of the k largest eigenvalues of a matrix M. Then the optimising 
V will be equal to the eigenvalues of the optimising p. To see this, recall that if v w, it holds 
i/(v) < Hi'w) [32| . Let us now assume that the (unique) optimising v of the Shannon entropy 
majorises but is not equal to the eigenvalues of the optimal p. Then Sip) would be greater or equal 
to Hiv), which is a contradiction. Since v w implies /(v) > /(w) for any convex function [32], 
the same argument can be applied to ( [ItI ). 

So far, we have replaced the von Neumann by the Shannon entropy and, in the dual problem, 
the matrix exponential by scalar exponential functions and transformed the arising eigenvalue 
constraints into semidefinite constraints. In order to apply Sedumi or SDPT3, all that is left to 
do is to approximate the objectives by piecewise linear functions. While this will only give us 
approximate solutions, let us note that the optimisers found in this way, can be easily checked to 
be feasible and inserted in the original objective. 

When using PPT or other positive maps the dimension of the matrix variables increases expo¬ 
nentially with the number of qubits. As we will show in the next section, this makes it difficult to 


go beyond five qubits. The computation of (16) greatly simplifies when only witnesses are used. In 


this case there will only be scalar variables p and Ui and no semidefinite constraints. This allows 
for application of a standard non-linear solver, such as FMINCON |33j . Using FMINCON we were 
able to obtain results for up to 13 qubits, as we will present in the next section. 

If we want to show if the Gibbs state is entangled, we will also need to compute the Gibbs state 
entropy Sipcibbs) = + InZ. To do so numerically it is sufficient to compute the eigenvalues 

of the Hamiltonian, which allows for computation of the partition function Z = Tr and the 

mean energy E = Ti pabbsH = 


D. Examples 

In order to test our method, we have implemented it for Heisenberg model, which was introduced 
in order to simplify the analysis of systems of spins, such as ferro- or antiferromagnets. It only 
takes into account the nearest neighbour exchange interaction between the spins as well an external 
magnetic field. In the one dimensional case, i.e. a chain of spins, the Heisenberg model is described 
by the following Hamiltonian 

N N 

^ = -Yl H ^ (Tf 

i=l i=l 


( 24 ) 




10 


where Jx, Jy, Jz are the coupling constants for the x, y, z-components of the spins, N the number of 
spins and B the external magnetic field in ^direction. a? denote the Pauli operators for the 

i-th spin. Let us assume periodic boundary conditions, i.e. a ring of spins. If Jx = Jy = Jz =• J, 
we are talking about an isotropic XXX Heisenberg model. J > 0 and J < 0 correspond to 
ferromagnetic and antiferromagnetic systems, respectively. If Jx = Jy ^ Jz, the system is called 
an XXZ system and so on. If only one component of the spin is considered, i.e. only one Jj ^ 0, 
the Heisenberg model reduces to the Ising model. 

The numerical results presented below have been obtained using either Sedumi or SDPT3 as 
well as Yalmip |34] . Using PPT constraints we have applied the method for up to five qubits, for 
witnesses for up to 13 qubits. 


1. Antiferromagnetic Heisenberg model 

For the one dimensional antiferromagnetic Heisenberg model (XXX with J = — 1), it is possible 
to detect entanglement at higher energies than UMl, both using the von Neumann and the linear 
entropy as proxies. The computations have been performed using both the partial transpose with 
respect to the partition Ai evenl^j odd and all possible partitions. 

Let us start with the results for the even versus odd partition. In tables|I]the energy ranges where 
entanglement is detected are shown for the von Neumann entropy and for the linear entropy, re¬ 
spectively. Here N denotes the number of qubits, Eq the ground state energy and F'min,PPT-even-odd 
denotes the smallest mean energy allowing for PPT w.r.t the even versus odd partition. Below 
that mean energy any state is guaranteed to be entangled. Note that Flmin.sep > ^'min,PPT-aii > 
-Emin,PPT-even-odd- Between Llmin,PPT-even-odd and Flmax,gap all states falling into the entropy gap 
are entangled. I.e. Flmax,gap is the largest energy for which the proxy method can work. This 
includes the Gibbs state in the von Neumann case. Note that the Gibbs state is not necessarily the 
maximiser for the linear entropy. Hence, instead of the Gibbs state we have computed the state 
with maximal linear entropy 


*S'L,max = max Slip) 

p state, ^Hp=E 


(25) 


and compared it to the constraint linear entropy. £ is defined as the fraction of the energy range 
where entanglement is detected 


_ .^maXjgap Eq 

F'max “ Eq 


(26) 


where Flmax denotes the largest energy eigenvalue. 

As can be seen it tables |IJ entanglement can be detected in a large area of the energy spectrum. 
The maximum mean energy where entanglement can be detected is also substantially higher than 
—N, which is the maximum mean energy where the methods of [Mnj work. Let us also note 
that the ground state energies given in [I 2 HII] are only correct in the limit of large N, as has been 
noted in m- The fact that the linear entropy detects more than the von Neumann entropy could 
be a result of the piecewise linear approximation of the von Neumann entropy in the maximisation, 
which results in weaker bounds on 5max,sep- 

The entropy gap can be seen in figure and[^ for the von Neumann entropy and for the linear 
entropy, respectively. Here the Gibbs state entropy or S'L,max) as well as Smax.pPT-even-dd are plotted 
versus the mean energy. Entanglement can also be detected in the presence of a magnetic field. 


Tables HI and IV show the detected energy ranges in a three qubit system for the von Neumann 
and linear entropies, respectively. The entropies are plotted in figure]^ for B = 3. 
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N 

Eo/N 

F^min,PPT-even-odd 

-E^max,gap / ^ 

£ 

3 

-1.000 

-1.000 

-0.610 

0.195 

4 

-2.000 

-1.000 

-0.660 

0.447 

5 

-1.494 

-1.008 

-0.695 

0.320 


TABLE 1. Energy ranges in which the von Neumann entropy can be successfully used a proxy for detecting 
entanglement in the antiferromagnetic Heisenberg model (J = —1). 


N 

Eo/N 

F^min,PPT-even-odd 

-E^max,gap / ^ 

£ 

3 

-1.000 

-1.000 

-0.600 

0.200 

4 

-2.000 

-1.000 

-0.370 

0.543 

5 

-1.494 

-1.008 

-0.302 

0.478 


TABLE 11. Energy ranges in which the linear entropy can be successfully used a proxy for detecting entan¬ 
glement in the antiferromagnetic Heisenberg model (J = —1). 




(a) E/N (b) E/N 

FIG. 3. The entropy gap for the antiferromagnetic Heisenberg model (J = —1) with 5 qubits and the PPT 
condition w.r.t. subsystems 1, 3, 5. (a) shows the von Neumann entropy, (b) the linear entropy. The Gibbs 
state entropy (red line) and the dual of the constraint, i.e. the witness canonical ensemble entropy (WCEE), 
(blue line) plotted versus the mean energy. The difference of the two lines constitutes the entropy gap. 


Going from PPT w.r.t. the even-uneven partition to all possible partitions greatly increases 
the computation time as the number of partitions grows exponentially with the number of 
qubit. Since the constrained entropy increases monotonically in E, it holds ii^max,gap,PPT-aU > 
-E'max,gap,PPT-even-odd) possibly increasing the energy range, where entanglement can be detected. 
In the example considered here, however, only a small increase can be obtained. See tables |V] 
as well as figure]^ For three qubits it holds F^ min ppT-aii/-/V = —0.6 > —1, implying that 


and VI 


-E'min,sep/-V > —0.6 > —1. This shows that the result of [l2l|Tl] is suboptimal for an odd number of 
qubits, as mentioned in M- The reason is that they use a partition into two sub-lattices such that 
every neighbouring sites belong to different sub-lattices, which is not possible for an odd number 
of qubits. 





























12 


B 

Eo/N 

-®min,PPT-even-odd/-^ 

^max,gap/ ^ 

£ 

0 

-1.000 

-1.000 

-0.610 

0.195 

1 

-1.333 

-1.333 

-0.720 

0.184 

2 

-1.667 

-1.667 

-1.033 

0.136 

3 

-2.000 

-2.000 

-1.540 

0.077 


TABLE III. Energy ranges for three qubits in the antiferromagnetic Heisenberg model (J = —1) in a 
magnetic field B = 3 using von Neumann entropy 


B 

Eo/N 

T'min,PPT-even-odd / N 

^max,gap / ^ 

£ 

0 

-1.000 

-1.000 

-0.600 

0.200 

1 

-1.333 

-1.333 

-0.560 

0.232 

2 

-1.667 

-1.667 

-0.717 

0.204 

3 

-2.000 

-2.000 

-0.940 

0.177 


TABLE IV. Energy ranges for three qubits in the antiferromagnetic Heisenberg model (J = —1) with a 
magnetic field B = 3 using linear entropy 


2. Dicke states 

Dicke states were first considered in the theory of coherent spontaneous light emission m- 
They are defined by 


\D^ ) = 



(27) 


where 



FIG. 4. The entropy gap for 3 qubits in the antiferromagnetic Heisenberg model (J = —I) with a magnetic 
field B = 3, again for PPT w.r.t. the even versus uneven partition, (a) shows the von Neumann entropy, 
(b) the linear entropy. 
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N 

Eo/N 

■£imin,PPT-all/^ 

-^max,gap / ^ 

£ 

3 

-1.000 

-0.600 

-0.600 

0.200 

4 

-2.000 

-1.000 

-0.660 

0.447 

5 

-1.494 

-0.809 

-0.809 

0.275 


TABLE V. Results for the von Neumann entropy in the antiferromagnetic Heisenberg model (J = —1) with 
a magnetic field B = 3 using all possible partitions 


N 

Eo/N 

-E’min,PPT-all/^ 

-£^max,gap / ^ 

£ 

3 

-1.000 

-0.600 

-0.594 

0.203 

4 

-2.000 

-1.000 

-0.360 

0.547 


TABLE VI. Results for the linear entropy in the antiferromagnetic Heisenberg model (J = —1) with a 
magnetic field B — 3 using all possible partitions 


m is the number of excitations, a denote sets of indices of excited subsystems and the sum is taken 
over all inequivalent sets of m indices. 

More recently Dicke states turned out to be a useful resource for quantum information processing 
task as they are LOCC transferable to GHZ or W states [36j. Several experiments have successfully 
created Dicke states, e.g. [36H38] . 

In [3^, it has been shown using perturbation theory that the approximate ground states of 
anisotropic ferromagnetic XXZ Heisenberg Hamiltonians (i.e. Jx = Jy ^ Jz > 0) are Dicke states. 
For B = 0 and > 0 the ground states are given by {|D^)}^^q with n + 1-fold degeneracy, 

which makes entanglement at low energies unlikely. This is because the system will soon go into a 
mixture of the Dicke states, which is separable. For Jx > Jz and B ^ 0, however, the degeneracy 
vanishes. In first order perturbation theory |D^) is the ground state for 


n — 2m + 1 
n — 1 


AJ <B < 


n — 2m — 1 


n — 1 


AJ 


(29) 


where AJ = Jx — Jz is the anisotropy parameter. 

A number of methods have been developed to detect genuinely multipartite entanglement in 



FIG. 5. The entanglement gap for all possible partitions for the von Neumann (a) and linear (b) entropies 
in the antiferromagnetic Heisenberg model (J = —1) with a magnetic field B = 3. 
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N 

m 

B 

AJ 

Eo/N 

-E^min,gap 

-E'max,gap / ^ 

T 

^ mm,gap 

T 

max,gap 

S 

11 

5 

-1 

10 

-13.753 

-13.599 

-12.890 

3.970 

7.930 

0.031 

11 

5 

-3 

20 

-26.707 

-26.299 

-24.698 

8.460 

16.420 

0.038 

11 

5 

-1 

20 

-26.525 

-26.344 

-24.802 

6.940 

15.850 

0.033 

13 

6 

-1 

12 

-16.276 

-16.149 

-15.577 

3.970 

7.930 

0.021 


TABLE VII. Detection energies and entropies for the XXZ 


Hamiltonian with = I using (30) as witness. 




(a) E/N (b) E/N 

FIG. 6 . Entropy gap plotted versus mean energy for (a) II qubits, m = 5, B = —3 and AJ = 20 and (b) 
13 qubits, m = 6, B = —1 and AJ = 12. 


Dicke states |4nH46) . of which |44j is the most generally applicable one. It consists of a non-linear 
witness, i.e. an non-linear inequality that has to hold for every biseprable state. In order apply the 
witness only a polynomial (in the number of qubits) number of local measurements is necessary. 

Combining [Sj with the use of proxy witnesses, entanglement can be detected even more easily 
experimentally. To do so, let us introduce a (weaker) linear version of the non-linear witness given 
in [Hj. It is given by 


IT’^ = - 

m 2 


^ ^ ( |d/3)(dQ,| |dQ,)((i^| -|- I dan/3)(^«n/31 T \daui3){daui3\) j 

(«,/3)e7 


(30) 


where 7 = {(a,/3) \af^ fj\ = m — 1}. detects \D'^) in a maximal way. Using these witnesses 
as a constraint in Q, it is possible to detect entanglement for different choices of Jx and Jz and 
B satisfying (29). Results for 11 and 13 qubits are presented in table VII and figures and 


III. TO THE THERMODYNAMIC LIMIT 
A. The ground state witness 

In order to extend the above results to the thermodynamic limit, let us now focus on witnesses, 
which are diagonal in the energy eigenbasis. A general example of such a witness is given as 
W = al — |Uo)(T^o|) where |iiio) is the ground state of H and a = maxj^^ \ {Eo\\(^)\ip)\‘^ [37]. 
Instead of |So)(-E'o|) one could also use any projector onto one or several energy eigenstates. The 
only obvious requirement being that the states are entangled themselves, otherwise the resulting 
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(a) i (b) A 

FIG. 7. Entropy gap plotted versus temperature for (a) 11 qubits, m = 5, B = —3 and AJ = 20 and (b) 13 
qubits, m = 6, B = —1 and AJ = 12. 


operators would not be witnesses. We can then either maximise over a fixed bipartition, or over all 
possible ones. In the former case it is possible to detect bipartite entanglement w.r.t. the partition 
chosen, in the latter case we can detect genuinely multipartite entanglement. Inserting the witness 
into (16), we see that the exponent becomes diagonal. 


v) = InTrexp {^Ei + u{a - 6io)) \Ei){Ei^ - ^E 
= In ^ exp {fiEi + v{a - 6io)) - ^lE, 


(31) 

(32) 


where Ei and \Ei) are the energy eigenvalues and eigenvectors, respectively and <5*0 is the Kronecker 
delta. Setting the gradient equal to zero, we obtain 

^£(^, u) = 0 =^{i) ^ exp {iiEi + v{a - (5io)) {Ei - E) = {) 

^ i 

d 

—£(^, I/) = 0 =>(fi) ^ exp {fiEi + u{a - 5io)) {a - 6io) = 0, 


which, by convexity of i, are sufficient conditions for a minimum. Solving those transcendent 
equations w.r.t. the Lagrange variables is only possible numerically for small N. Still it is possible 
to obtain a result for the thermodynamic limit, as we will now demonstrate. 


Let us first consider the case where a = 1, i.e. H has a separable ground state. Since the first 
term in the sum of (ii) vanishes, (ii) cannot be fulfilled for any > 0. Since i increases with v, the 
minimum is attained at = 0. At = 0, however, (16) is an upper bound on the unconstrained 
problem, namely 


minln y^exp(^£'i) - ^lE > S{pGibbs) 

I 


Since S{pGibbs) = lnZ + I3E^ the minimum is attained sX p, = —/3. In this case no entanglement is 
detected. 

Let us now move on to the case where a < ^ ^ ° . Assuming that the minimum is attained at 
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u = 0, i.e. that (ii) is fulfilled at u <0, (ii) tells us that 


exjp{—/3Eo + v{a — 1)) _ e 

Xli exp(-/3£lj + v{a - Jio)) Y.i 


g— /3Eo g— /3Eo 

g-^Eo _|_ “ Z 


with equality for i/ = 0, which is a contradiction. Hence the minimum is attained at > 0. This 
implies that min,,>o,^^< fi) = S{pGibbs), hence 5'max,AA,lVi,E < S{pGibbs)- 


Theorem 4 Let SEP be the set of separable pure states w.r.t. some partitions anda = max^fp'f^sEp\{Eo\\4>)\‘^ ■ 

Then, if a < ^ , the Gibbs state as well as any state with entropy Smax,AA,Wi,E < S < S{pGibbs) 

will be inseparable w.r.t. those partitions. 

Note that while for a < ^ ^ ° entanglement of the Gibbs state can also be detected by applying the 
witness directly, the same is not true for states with entropy *S'max,AA,Wi,£; < S < S{pGibbs)- Since 
such states always exist, this shows that our method generically works. Theorem]^ can be applied 
in the thermodynamic limit. Since Eq generally scales linear in N while Z does so exponentially. 
Theorem]^ can yield non-trivial results where a = 0{e~^). 


B. Example 


We consider a special case of the Heisenberg Hamiltonian, given by 


N 




i=l 


1 + r 


I 

^i+i + 


1 — r 


+ hat 


This is known as the XY model in a transverse magnetic field. If r = 1, we talk about the Ising 
model, if r = 0, the XX model. The XY model undergoes a phase transition at /i = 1 [2]. The 
ground state energy as well as the partition function have been computed in the thermodynamic 
limit by [IH]- According to a conjecture numerically tested in [3^, a w.r.t. full separability is given 
by 


In a , 1 

hm -= 2 max / auln 

N^oc N i Jo 


where tan20 = f[_JJtys 2 nii — f < 0 < f • Using this, we can numerically show entanglement for 
a wide range of T and h. See Figure]^ In particular we can show entanglement in vicinity of the 
phase transition, in accordance with [2]. 


C. Arbitrarily large Gap 

Unfortunately the entropy gap closes in the limit when using the regular ground state witnesses, 
a fact which is related to the lack of robustness of the ground state witness itself. It would be 
most interesting to find other, asymptotically robust witnesses, which could result in a macroscopic 
entropy gap. Let us now present an artificial example of a Hamiltonian which, with the right choice 
of witness, allows for an arbitrarily big entropy gap in the thermodynamic limit, thus proving the 
feasibility of this endeavour in principle. 
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Ising 



FIG. 8. Our results for the (a) Ising, (b) XY with r = 0.5 and (c) XX models. The l.h.s. show the maximum 
temperature at which Theorem yields non-trivial results. The r.h.s is a contour plot of the positive part 
of limjv_>.oo ~ J- We can show entanglement where this expression is positive, i.e. left of 

the zero line. 
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The Hamiltonian is given by 

d^-i 

H =Y,k\^k){^k\, (33) 

k=0 

where d = 22 and are a Bell state basis of an n-qubit Hilbert space. The basis is 

ordered in such a way that YliZo + for ^ = 0...d— 1. Note that these first d basis 

elements have orthogonal support. As witness, we choose 


d-l 


W = al-Y,\'^k){-^k\, 


(34) 


k=0 


which, again, is diagonal in the energy eigenbasis. a is the defined as the maximal overlap of the 
projector a separable state. As all the |'hfc) have orthogonal support, a is 

the maximal overlap of any of the IT^) with a separable state, hence a = ^ because the |Tfc) are 


maximally entangled. Inserting this into (19), we obtain 


AS > In 


Y^d^-l pk 

l^k=0^ 




(35) 


ed 


Letting n, hence d, go to infinity, we obtain 


lim AS > u, (36) 

n^oo 

where we have used the fact that lim^^oo Ylk=o^ = lim(^_^oo Ylk=o Note that v 

can be chosen arbitrarily large. 


IV. DISCUSSION 

A. Robustness 


In macroscopic systems assuming the exact form of the Hamiltonian is always an idealisation. 
While it is fair from a physical point of view one might wonder about the impact of mis-characterised 
Hamiltonians on onr entropy gap AS. If we assume that the real Hamiltonian is given as H = 
H + eP we can that the dnal entropy of the actnal Hamiltonian for the chosen parameters 

<Smax,in,E(Ai, k') = log Tr exp{-fj.H - p.eP + vW) + pE + ^eTr(pP), (37) 

can be bounded from above using the Golden-Thompson Tr(e^+^) < Tr(e"^)Tr(e^) and Holder’s 
||Ai?||i < ||A||oo||B||i inequalities to yield 

Sraa.^,W,E{fJ-, e) < 5'max,W,£ (li, e) + ;UeTr(pP) - pePo ■ (38) 

This shows that small perturbations or inaccuracies in the description of the Hamiltonian will only 
have a correspondingly small impact on the validity of the entanglement certification by proxy 
witnesses. 
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B. Experimental estimation of entropy 

While the mean energy and the entropy are both macroscopic properties of quantum states, 
measurements of the latter are not possible directly as they do not correspond to a quantum 
observable. There are however various ways that the global entropy can be determined, and all 
we need is a lower bound on the entropy. The most straightforward way would of course consist 
of equilibrating the system with a thermal bath at temperature T. This is naturally the case, as 
any system found in nature at ambient temperature T, will most accurately be described by its 
corresponding Gibbs state. Thus, once the system is equilibrated we know that its corresponding 
entropy should correspond to the Gibbs entropy S. Starting from equilibrated systems 
one can introduce global quenches of system parameters, such that the entropy will still be bounded 
from below by the initial Gibbs entropy due to the second law of thermodynamics, despite the 
system being far out of equilibrium. In that way one directly receives a lower bound on the system 
entropy for a wide range of out of equilibrium systems and experimental preparations. There 
are of course other methods, such as reasonable assumptions about symmetry in the state, that 
can be used to estimate linear entropy from macroscopic spin observables (as would be possible 
e.g. in Ref. [50]). Another straightforward way to obtain lower bounds would be access to a 
d X d sub-matrix of the global state, whose entropy will always yield a lower bound to the global 
entropy. The exact ways of experimentally estimating entropies is of course highly dependent 
on the experimental setup, access and reasonable assumptions about system properties. As any 
non-zero amount of entropy already provides an advantage in entanglement detection it is however 
fairly straightforward to infer correspondingly useful entropies. 

This highlights the main advantage of our approach: While entanglement of Gibbs states could 
be directly inferred from its description, it would strictly work only if we know that the state is in 
fact in thermal equilibrium. This would require strict assumptions about the state of the system 
and exact characterization of its Hamiltonian. Using our entropy gap there is no need whatsoever 
to assume any particular form of the system’s state: measurements of the mean energy and lower 
bounds on the global entropy are completely sufficient to prove that the underlying system is 
entangled in a robust way. We believe that future work will uncover witnesses robust enough to 
be amenable to the proxy method, yielding a sizeable entropy gap in the thermodynamic limit. 


V. CONCLUSION 

We have introduced a framework in which concave functions can be used as a proxy for detecting 
entanglement in many body systems. As the most relevant example we have explored the use 
of entropy in this context, yielding separable Gibbs ensembles and a corresponding entropy gap 
that can be used for entanglement detection. Through the efficient description through witness 
canonical ensembles this unlocks powerful tools from entanglement theory, such as positive maps or 
entanglement witnesses to be harnessed in situations where the actual estimation of said quantities 
is experimentally impossible. 

We hope that our work can contribute to the further understanding of the relationship between 
important physical aspects of many-body systems, such as phase transitions, and the paradigmatic 
feature of Quantum Information Theory, entanglement |3]. We address some possible new avenues 
in entanglement detection, by indicating the importance of an open conjecture regarding the geo¬ 
metric measure of entanglement for ground states |l9] and by specifying particularly useful forms 
of entanglement witnesses in many-body systems. 

We have furthermore demonstrated the usefulness and feasibility of this approach in exemplary 
and paradigmatic physical Hamiltonians. For moderate system size we have demonstrated how 
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our results directly improve upon previous work on inferring entanglement from macroscopic ob¬ 
servables. The resulting entropy gaps quantify to what extent even out-of-equilibrium systems can 
be certified to exhibit entanglement. In the thermodynamic limit we have shown that particularly 
promising paths towards this goal are strongly connected to ground state properties of many-body 
Hamiltonians. Our results in this context furthermore elucidate how a development of robust en¬ 
tanglement witness techniques can prove useful, even if they themselves remain inaccessible due to 
experimental limitations. 
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